Field-induced Ordering in Critical Antiferromagnets 
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Transfer-matrix scaling methods have been used to study critical properties of field-induced phase 
transitions of two distinct two-dimensional antiferromagnets with discrete-symmetry order parame- 
ters: triangular-lattice Ising systems (TIAF) and the square-lattice three-state Potts model (SPAF- 
3). Our main findings are summarised as follows. For TIAF, we have shown that the critical line 
leaves the zero-temperature, zero-field fixed point at a finite angle. Our best estimate of the slope 
at the origin is (dT c /dH) T=H=0 = 4.74 ±0.15. For SPAF-3 we provided evidence that the zero-field 
correlation length diverges as £(T — > 0, H — 0) ~ exp(a/T x ), with x = 1.08 ±0.13, through analysis 
of the critical curve at H 7^ plus crossover arguments. For SPAF-3 we have also ascertained that 
the conformal anomaly and decay-of-correlations exponent behave as: (a) H = 0: c = 1, 77 = 1/3; 
(b) H ± 0: c= 1/2, t) = 1/4. 

PACS numbers: 05.50.±q, 05.70.Jk, 64.60.Fr, 75.10.Nr 
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I. INTRODUCTION 

Frustrated systems with macroscopically degenerate 
ground states still pose intriguing theoretical challenges. 
For pure Ising and g-state Potts antiferromagnets (AF's), 
in particular, this degeneracy results solely from the in- 
terplay between "dynamics" (i.e., the number of states 
per lattice site) and geometry (lattice topology). In this 
paper we shall deal exclusively with Ising and 3-state 
Potts AF's, respectively on the triangular and square 
lattices. It is well established that, in the absence of 
an external magnetic field, both systems display a criti- 
cal ground state (in the sense that spin-spin correlations 
decay algebraically with distance) , and are paramagnetic 
for all temperatures T > |l]-|^]. In both cases, a uni- 
form field H removes the residual entropy per spin, and 
long range order can set in at finite temperatures, below 
a field-dependent phase boundary T C (H) 

For the Ising antiferromagnet on the triangular lattice 
(TIAF), the phase transition at T C {H) has been deter- 
mined to be in the 3-state ferromagnetic Potts model 
universality class ^,^] ; subsequent analysis in the context 
of conformal invariance led to a conformal anomaly 
(or central charge) c = 4/5, which is also consistent 
with the value for the 3-state Potts ferromagnet [fl0|l . 
Further, it has been found that the critical phase at 
T = H = extends into a small region T = 0, H < 
H KT ^ 0.27 that is, H is not a relevant scal- 

ing field at T = 0, as initially thought ||. At Hkt 
the system undergoes a Kosterlitz-Thouless (KT) transi- 
tion to a long-range ordered state p3| . Within the zero- 
temperature critical phase, one has continuously-varying 
critical exponents and, accordingly the conformal 
anomaly is c = 1 |l3| . Much less is known about field ef- 
fects on the 3-state Potts antiferromagnet on the square 
lattice (SPAF-3), apart from indications that, for H ^ 
the transition at the corresponding T C (H) belongs to the 
two-dimensional Ising model universality class 



Since the scaling behaviour (i.e., relevance or marginal- 
ity) of the uniform field at T = influences the shape of 
the critical curve near T = H = in a fundamental 
way, an accurate evaluation of T C (H) is clearly of inter- 
est. With this in mind, here we investigate the finite- 
temperature field-induced transition in both the TIAF 
and the SPAF-3, by means of transfer-matrix scaling 
methods |l5],|l6| • Fot TIAF we concentrate on the shape 
of the critical curve close to T = H = 0, for reasons 
to be stated in the corresponding section. For SPAF-3 
we determine the critical curve, as well as the conformal 
anomaly and the decay-of-correlations exponent 77 along 
it. The layout of the paper is as follows. In Sec. || wc 
outline our calculational procedure for the free energy 
and the correlation length, from which we determine the 
conformal anomaly and the critical curve, respectively. 
Results for the TIAF and SPAF-3 are presented in Sec- 
tions [II and IV, respectively. 



findings. 



Sec. ^ 



summarizes our 



II. MODELS AND TRANSFER MATRIX 
SCALING 

We consider infinitely long strips of width L, with pe- 
riodic boundary conditions in both directions. Ising or 
Potts spins sit on lattice sites and interact with each 
other, as well as with a uniform field, according to the 
Hamiltonian (including the multiplicative factor — /3 = 
-1/AfeT), 



(1) 



(id) 



where the first sum runs over nearest neighbour sites of 
either a triangular or a square lattice, depending, respec- 
tively, on whether <7j is taken to be or 1 (Ising), or 0,1, 
or 2 (3-state Potts); a convenient strip geometry for a 
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triangular lattice corresponds to the usual square strip 
with additional bonds along a fixed diagonal direction. 
K is the exchange coupling constant and the field H has 
been taken along the O-direction. 

We shall use units in which, for the triangular Ising 
/erromagnet, T" 1 = K c = ±ln-\/3, and for TIAF the 
upper critical field [such that T C {H > H c ) = 0] is H c = 
6 (6|9J. For SPAF-3, the corresponding quantities are 
T- y = ln(V3+ 1) (ferromagnet), H c = 4 §. 

As usual |l5[|l6t| , the free energy per spin /l (T, H) and 
the correlation length £l (T, H ) are given by 



h(T,H) 



InAi; H\T,S) 



-Cln(A 2 /A!) (2) 



where Ai (A 2 ) is the largest (second-largest) eigenvalue 
of the transfer matrix between two successive columns; 
the geometric factor ( = 2/V3 for triangular, and 1 for 
square lattices. 

We have obtained finite-size estimates of the criti- 
cal line by standard phcnomenological renormalization 
(PRG) procedures flD||l6| |: for fixed H we consider pairs 
of strips of respective widths L and L' and solve the im- 
plicit equation 



LZZ 1 {T\H) = L' &(T*,H) 



(3) 



for the fixed-point temperature T*. This approach is 
rather safe because the only underlying assumption is 
that a second-order phase transition occurs, without any 
further hypothesis on its universality class. As explained 
in detail below, we will be particulary concerned with the 
proximity of (T = 0, H = 0), where crossover between 
different sorts of critical behaviour is expected. Owing to 
sublattice symmetries, L and L' must be both multiples 
of 3 for TIAF (where we use V = L - 3, L = 6, . . . , 18), 
and 2 for SPAF-3 (respectively L' = L-2, L = 4, . . . ,12). 
At some special points, such as (T = 0, H = 0), we went 
up to L = 14 for SPAF-3. 

For each fixed H the sequence of finite- L estimates of 
T* was extrapolated against suitable inverse powers of 
L, so that a set of temperatures T^ xtr (H) was produced, 
which represents our best approximation to the true crit- 
ical curve. We then calculated the free energy and the 
correlation length for finite L, as in Eq. (||), at T^ xtr {H). 
From these we produced finite-size estimates of the con- 
formal anomaly c and the decay-of-correlations exponent 
77, respectively via 



L 2 [f L {T c 



and 



foo{T c )] — g , 



Ltt 



V 



(4) 



(5) 



as given by conformal invariance [T^Jij], where the in- 
dependence of T c (and thus, of c and 77) is implicitly un- 
derstood. The sequences of finite- L estimates cl and tjl 



were again extrapolated to L — > 00 to give the final values 
c(H) and r)(H). 

The extrapolated phase boundaries were usually ob- 
tained under the simplifying assumption of single-power 
corrections to scaling: for given H we assumed 



T*(L,H) -T C (H) - IT* 



(0) 



It is expected [ pL 8[Jig[ | that ip = oj + X/v, where u> is the 
leading correction-to-scaling exponent: 



£ L (T c )=AoL (l + A\L~ 



(7) 



and v is the correlation- length critical exponent, known 
to be 5/6 (1) for 3-state Potts (Ising) ferromagnets [rel- 
evant for TIAF (SPAF-3) in non-zero field]. However, u> 
is not known in advance: though in several cases pc| ] nu- 
merical evidence has been given in support of to — 2, this 
is in principle a non-universal quantity which depends on 
the pertinent operator algebra, and probably on lattice 
effects as well; e.g., for Ising ferromagnets on triangular 
and honeycomb lattices, it has been found |2l]] that uj = 4 
fits Eq. (M) extremely well. Attempts to keep u> fixed at 
2 in Eq. (n) for the present case resulted in fits of widely 
varying quality. Thus we took a pragmatic view, and for 
each H varied ip within reasonable limits (to be spelt out 
below) until a good fit turned out. More often than not, 
the smallest-L term [L = 6(4) for TIAF (SPAF-3)] was 
discarded. Typical uncertainties for T^ xtr (H) were one 
part in 10 4 , which means that the dominant contributions 
to final spreads in i](H) and c(H) are attributable to the 
respective extrapolations of r]L and Cl to L — > 00 . From 
Eq. (0) , the corrections to r)L of Eq. (||) are expected to 
behave as . Surprisingly, a fixed uj = 2 gave reason- 
ably good fits throughout the range of fields investigated, 
for both models. For the conformal anomaly the addi- 
tional unknown foo(T c ) arises; assuming corrections to 
scaling to Eq. (Q) also with w = 2, (in this case, such 
corrections have been shown to work well for Potts p2| 
and Ising [^3|] ferromagnets) we performed least-squares 
fits of our data to a parabolic form in L~ 2 p3]. 

Further, the field dependence of c and 77 can be anal- 
ysed within a finite-size scaling (FSS) theory of crossover 
effects p3. We first assume the existence of two bulk 
correlation lengths, £° and £ T : the former is only di- 
vergent at T — > 0, H = 0; the latter diverges both at 
T — > 0, H = (i.e., £° ~ £ T in this case), as well as 
at T C (H) > 0, with different asymptotic forms. The two 
scaling variables are then and L/£ T , which allows us 
to cast the finite-size correlation length and free energy 
in the forms 



and 



Hl(T,H) = LQ(L/e,L/e), 



Il(T, H) = L~ d R{L/e,L/?), 



(8) 



(9) 



respectively, where Q and R are extended scaling func- 
tions. The crucial difference between the scaling of these 
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and of any other quantity (e.g. susceptibility, specific 
heat, magnetization) is that the leading power in the In- 
dependence is fixed, instead of a ratio of critical expo- 
nents, x/v, which may change from xq/vq to xt/^t- In- 
deed, the main IZ-dependence in the crossover function is 
expected to arise as [f(H)] e , with e = xt/vt~ xq/vq |p4) . 
Thus, low-order corrections to the asymptotic behaviour 
must be wiped out; the field-dependent crossovers in both 
77 and c are therefore expected to be quite fast. We shall 
see that these predictions are borne out rather well by 
numerical data for SPAF-3. For T1AF, technical difficul- 
ties (to be described) connected to extrapolation of the 
critical boundary translate into a more mixed picture. 



III. TRIANGULAR ISING ANTIFERROMAGNET 

For TIAF in the range 1 < H < 5.5, best fits to Eq. 
(|J) were attained with ip ~ 3.5 — 5.5 (higher values for 
lower fields). In that region our PRG estimates extrap- 
olate to values virtually identical to those found in Ref. 
||. Those authors started from the assumption that, for 
all H ^ 0, the TIAF is in the same universality class 
as the 3-state two-dimensional Potts ferromagnet MM, 
and located the points (T, H) where £c(T, H) = IbL/Air, 
corresponding to 7/ = 4/15 as given by conformal invari- 
ance p7| . Procedures of this sort were put forward by 
Blote and den Nijs [Q, and are expected to be less vul- 
nerable than PRG to numerical inaccuracies, provided 
the universality class of the transition is not in doubt. 
However, near a multicritical point (such as T = H = 
here) crossover effects may also take their toll. Indeed, 
even though there is no a priori reason to question uni- 
versality in the present case, convergence of the data of 
Ref. H deteriorates rapidly close to the origin, to such an 
extent that the authors quote no extrapolations for the 
critical line for H < 0.5. In our investigation, we found 
that for 0.15 < H < 0.5 the PRG curves crossed each 
other, thus making the extrapolation procedure unwork- 
able. An example can be seen near the right edge of Fig. 
HI below, which also shows that closer to the origin the 
curves again behave monotonically against L. 

Before giving details of extrapolation in that region, we 
recall that the shape of the low-T, low- if phase boundary 
was discussed in Ref. || . At the time it was believed that 
H was a relevant scaling field along T — 0, from which it 
was concluded that, since the correlation length diverges 
as exp(a/T) |2j| at H = 0, T -> 0, the critical line should 
approach the origin tangentially to the T- axis. On the 
other hand, PRG with L = 6 and 9, in the notation of 
our Eq. (||) , yielded a finite slope at the origin. This was 
interpreted as a deficiency of the calculational method || . 
Here we have reexamined the matter by extending PRG 
to L = 18, extrapolating our data and making contact 
with the more recent results |[2],|l3| which point to the 
existence of a Kosterlitz-Thouless phase at T = 0, H 7^ 0. 
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FIG. 1. PRG estimates of the critical line of TIAF near the 
origin, obtained by solving Eq. L values given on figure 
(L' = L-3). 



In Figure [j] the fixed-point solutions of Eq. (|J) are dis- 
played. Though convergence becomes prohibitively slow 
for H < 8 x 10 -3 , there is plenty of leeway to establish 
that, as H — > 0, all our finite-!/ curves become straight 
lines which (to within one part in 10 4 — 10 5 ) cross the 
origin. The straight sections become shorter with grow- 
ing L, thus one must be careful before predicting a finite 
slope at the origin for the actual phase diagram. In Table 
I we show the slopes Sl of the straight sections of finite- 
ly curves, as well as the values H max (L) above which 
said curves begin to deviate from linear behaviour. Error 
bars for H max (L) are somewhat subjective, but certainly 
quite conservative, as can be seen from visual inspection. 
Should L — > 00 extrapolation produce a definitely neg- 
ative value of H max , one could be sure that the finite 
slope is a finite-size artifact. However, we have found 
H max (oo) = 0.01 ± 0.02 from a rather good scaling of 
our data against L _1 . 

Bearing in mind that the only "small" typical field nat- 
urally arising in the problem, Hkt, is one order of mag- 
nitude larger than this (thus a strictly positive H max (oo) 
of order 10~ 2 would have no clear physical origin), we 
interpret the above result as signalling that H max (pd) is 
exactly zero. So, (i) the critical line does leave the origin 
at a finite slope (for which our best estimate, 4.74±0.15, 
comes from extrapolation of the Sl against L~ 2 ); but 
(ii) the straight-line part of the critical curve is of zero 
extent: one only has d 2 T c /dH 2 — at the origin. This 
latter quantity must be negative for all H > 0, as no 
inflection points are expected. 

These conclusions are consistent with the presence of a 
critical phase on the H- axis near the origin. The expo- 
nentially diverging correlation length at H — 0, T — > 
is roughly in balance with the (already infinite) ^kt 
along the zero-temperature axis. This way the critical 
line, where crossover between temperature- and field- 
dominated behaviour takes place, starts at finite angles 
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with both axes. 

Going back to extrapolation of the critical line for 
< H < 0.15, we first note that, by construction, 
our procedure of fixed-i? extrapolation automatically 
yields a straight line with slope 4.74 ± 0.15 for all H < 
H max (18) ~ 0.09. From the preceding arguments on the 
extent of the straight-line part, and concavity, of the crit- 
ical curve, this is an upper limit: T^ eal (H) < T^ xtr (H). 
Also, for 0.09 < H < 0.15 we have not managed to pro- 
duce good fits with single powers; instead, we were forced 
to resort to two-power fits using L^ 1 and L~ 2 . These 
facts have strong effects on the evaluation of c and r\ 
near H = 0, which we now turn to discuss. 
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FIG. 2. Conformal anomaly c (upper curve) and exponent 
r\ (lower curve) along the extrapolated critical line of TIAF 
near the origin. Expected values are (a) H — 0: c = 1, 
7] = 1/2; (b) H ^ 0: c = 4/5, -q = 4/15. 



Recall that at T = H = one has |25j,|lJ c = 1 
and r] = 1/2, while for H ^ the 3-state Potts values 
c = 4/5, r\ = 4/15 are believed to hold B. Our results 
for c and r\ along the extrapolated critical line, near the 
origin, are shown in Figure @. From the discussion in 
Sec. [n] one might assume that, apart from higher-order 
crossover effects, both quantities should behave in a step- 
function fashion. On the contrary, we see that they hover 
around their zero-field values for a significant range of H, 
which coincides with that where our extrapolation gives 
a straight line. Further on along the H- axis, conver- 
gence begins to deteriorate. Taking into account (i) the 
conclusion that straight-line sections of PRG curves are 
finite-size effects; (ii) the inescapable distortion imposed 
by them onto our constant--^ extrapolations; plus (iii) 
the fact that the true critical line is only expected to be 
straight at the origin, where it is joined by the KT line, 
we tentatively interpret the plateau-like behaviour of c 
and i) as a manifestation of the KT phase in an artificial, 
finite-size-induced fashion. We have not yet managed to 
propose a numerical test of this idea; however, as shown 
in the next Section, a measure of self-consistency of the 



argument is found in SPAF-3, where both the KT phase 
and the anomalous behaviour of c and r\ are absent. Note 
also that, on general grounds, estimates of L/tt^(T, H) 
at T < T^ xtr (H) (as T r c eal must be) would certainly pro- 
duce numerical values smaller than those displayed as 
rj in the Figure, which is not inconsistent with the ex- 
pected T) — 4/15. Finally, as regards 0.10 < H < 1.0 
(at which upper extremity one already has the results of 
Ref. H), the above-mentioned difficulties with L — > oo 
extrapolations of PRG curves translate into unsurmount- 
able obstacles to estimations of c and rj. 



IV. SQUARE LATTICE 3-STATE POTTS 
ANTIFERROMAGNET 

We begin the discussion of SPAF-3 by examining the 
point T = H ~ 0, where strips of maximum width L = 14 
sites were used. We have calculated c and 77 and found, 
after extrapolation, 

c = 0.999 ±0.001 ,77 = 0.333 ±0.001 ,(T = H = Q). 

(10) 

Owing to the unusually slow convergence of finite- L data 
for c in this case, we formed three-point fits with the 
sets {fi(T c )}, I — L, L — 2, L — 4. The sequences of ex, 
each estimate resulting from a three-point fit, were then 
extrapolated by a Bulirsch-Stoer algorithm f2(|^7j , which 
essentially amounts to assuming a single-power correction 
to scaling. Our best fit corresponded to that power being 
around 2. 

The unitary value of c has been predicted for the 
case [ p8p9[ | and is consistent with continuously-varying 
critical exponents fll4|], thus one might expect e.g. a KT 
phase on the T = axis, by analogy with TIAF. We 
shall return to this point below. Our result for rj ap- 
parently contradicts the direct evaluation of correlation 
functions of Ref. |oj, which yields rj = 1.33±0.02. To ex- 
plain this, we recall the prediction of Ref. Q which, for 
SPAF-3 with only first-neighbour interactions (in their 
language: v = 1, fi = 2n/3, yx = 1/2), reads: 



G(r) 



A B 

r 4/3 r l/3 



(11) 



where G(r) is the critical (T = 0) spin-spin correlation 
function at distance r; the sign of the second term de- 
pends on whether the two sites are on the same sublat- 
tice Q, that is, it is associated to the staggered magneti- 
zation. The result of Ref. |}0| is for correlations between 
spins on opposite corners of N x N finite lattices with 
N < 15, thus it gives the decay of the uniform magne- 
tization, dominated at short distances by the first term 
of Eq. ([ll]). Indeed, a plot of their data in the form 
r 4 / 3 G(r) against r appears to approach a straight line 
for large r, with A ~ 1.4 and B ~ 5 x 10~ 3 . On the 
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other hand, by relying on the amplitude-exponent rela- 
tion given by conformal invariance (l7| , our approach au- 
tomatically picks up the behaviour of the smallest gap (or 
longest correlation length) of the transfer matrix, which 
indeed couples to the staggered magnetization. These 
considerations were very recently rederived via a height 
representation of the model [ pl| and confirmed by nu- 
merical work |32|| . In the present case, our estimate is 
entirely consistent with the second term of Eq. ([ll]), and 
also with Monte Carlo work |32|;[33| . 

We have paid special attention to the shape of the 
critical curve near the origin. Throughout the range 
0.002 < H < 0.011, we managed good fits to Eq. (|) 
with ip in the range 4.9 — 6.7. For now, we concentrate 
on the analysis of that region. 

We recall that, although it is agreed that in zero field 
the system is critical only at T = 0, there seems to be 
no consensus about how the correlation length diverges, 
except in that an exponential singularity £(T — > 0,H = 
0) ~ exp(a/T x ) is present. The value of x has been vari- 
ously estimated as 1.3 (by analysis of the Roomany-Wyld 
approximant 34 in a transfer-matrix calculation 0, and 
Monte Carlo (MC) work Q); 1 (further Monte-Carlo 
work p5|) and 3/4 (conformal invariance arguments cou- 
pled with an analysis of the eigenvalue spectrum of the 
transfer matrix |28[). 

If we assume that, along the T = axis, H is a rele- 
vant variable with scaling index yn , a standard crossover 
argument j(| implies that on the critical curve T c ~ 
| InH]^ 1 ^. If, on the other hand, an extended criti- 
cal phase is present as in TIAF, the results of Section 
[II indicate that such shape is unlikely to be found. To 
be fair, we must point out that there is no compelling 
symmetry-based argument (such as vortex unpinning for 
TIAF [|l2|,|l3|) which leads one to infer the possible exis- 
tence of a soft phase here. 

In Figure || we show our results for T c \ \a.H\ x / x against 
H for x = 3/4, 1, 4/3 as well as our best fit for an asymp- 
totically horizontal line as H — > 0, which corresponds to 
x = 1.08 ± 0.13. This estimate and its respective error 
bar are based on analysis of the insert of the Figure: the 
x = 1.21 curve flattens at H ~ 0.002 (the lowest field we 
can reach), so those for x > 1.21 certainly bend down- 
wards before touching the vertical axis. Analogously, the 
x = 0.95 curve is roughly straight, so those with x < 0.95 
will be concave upwards. The central estimate is taken 
as the average of these upper and lower limits. 

Thus we conclude that (i) there is no numerical evi- 
dence of an extended critical phase at (T — 0, H ^ 0); 
(ii) our data are consistent with an infinite slope of the 
critical curve at the origin, meaning (via the crossover ar- 
gument above) that (in) the zero-field correlation length 
diverges with T -> as £(T -> 0, H = 0) ~ exp{a/T x ), 
x = 1.08 ±0.13. 

Of all previously available estimates, the latter value 
of x is only consistent with the Monte Carlo results of 
Ref. p5| . Those authors mention possible logarithmic 
corrections, which would give an enhanced effective ex- 



ponent. This would also be in line with the fact that our 
central estimate is slightly above unity. 
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FIG. 3. Plots of T c \ lnH\ 1/x for SPAF-3 for x = 3/4, 1, 4/3 
and 1.08. Insert: behaviour of plots for x = 0.95, 1.08 and 1.21 
from which our central estimate and respective error bars, 
x — 1.08 ± 0.13 have been extracted (see text). All curves 
normalized to one at H = 0.010. Here T c stands for (L — > oo) 
extrapolated values. Error bars coming from extrapolation 
are smaller than symbol sizes. 



Further on along the ii-axis, for 0.013 < H < 0.18 
the PRG curves for L = 8, 10, 12 crossed each other at 
nearly zero angle. In that region, we simply took straight- 
line fits of the three respectively values of T* against 
to obtain T^ xtr . However, for H > 0.2 monotonic 
behaviour returned, once again allowing use of Eq. (j^) 
with if> ~ 2.1 — 5.9. Additionally, along either sections 
of the extrapolated critical line matched one another so 
well across the gap, that they are joined by continuous 
lines 

In Figure |i| estimates of both c and 77 along the ex- 
trapolated critical line are displayed. One can see that 
for both quantities, the somewhat ad hoc extrapolation 
procedure in the intermediate- H region produces sensi- 
ble estimates, which join the adjacent sequences rather 
smoothly. 

Further, this time the predictions of Section || are seen 
to hold: apart from higher-order crossover effects, both 
quantities behave close to step-functions, converging to 
the respective Ising values c = 1/2, 77 = 1/4. At H = 
0.5 the curve for c has a minimum. We used strips of 
width L < 14 to produce an accurate estimate both of 
T c and c, which turned out as c = 0.47 ± 0.002. Thus we 
conclude that the (unaccounted for) residual crossover 
effects produce deviations of order at most ~ 6%. 
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FIG. 4. Conformal anomaly c (upper curve) and exponent 
77 (lower curve) along the extrapolated critical line of SPAF-3. 
Expected values are (a) H — 0: c = 1, 77 = 1/3 (see text); (b) 
H=£0: c=l/2, 77 = 1/4. 



V. CONCLUSIONS 

We have studied critical properties of field-induced 
phase transitions of selected two-dimensional antifer- 
romagnets with discrete-symmetry order parameters. 
Throughout our work, we attempted to minimise numer- 
ical effects originating from crossover between different 
universality classes, by applying carefully selected proce- 
dures both for finite-size calculations and for extrapola- 
tion of finite-size data to the infinite-lattice limit. For 
TIAF we did not entirely succeed, owing mainly to resid- 
ual effects ascribed to a Kosterlitz-Thouless phase along 
the zero-temperature axis. For SPAF-3, where our evi- 
dence shows that such phase is not present, we present 
results which are clean and unambiguous, for all quanti- 
ties investigated. 

Our main findings are summarised as follows. For 
TIAF, we have shown that the critical line leaves the 
zero-temperature, zero-field fixed point at a finite an- 
gle. Our best estimate of the slope at the origin is 
{dT c /dH) T=H=a = 4.74 ± 0.15. For SPAF-3 we pro- 
vided evidence that the zero-field correlation length di- 
verges as £(T -> 0,H = 0) ~ exp(a/T x ), with x = 
1.08 ± 0.13, through analysis of the critical curve at 
H 7^ plus crossover arguments. For SPAF-3 we have 
also ascertained that the conformal anomaly and decay- 
of-correlations exponent behave as: (a) H — 0: c = 1, 
77 = 1/3; (b) H ^ 0: c = 1/2, 77 = 1/4. 
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TABLE I. Slopes Sl and upper limits H max (L) of 
straight-line portions of approximate (PRG) critical lines for 
TIAF. From data displayed in Figure [j]. Extr. stands for 
extrapolated as L — > oo (see text). H max {6) is omitted, as 
extrapolation only took L — 9 — 18 into account 



L 


Sl 




6 


1.4979 ± 0.0001 




9 


2.0414 ± 0.0002 


0.170 ± 0.010 


12 


2.4564 ± 0.0002 


0.130 ± 0.010 


15 


2.8078 ± 0.0002 


0.105 ± 0.010 


18 


3.1200 ± 0.0010 


0.090 ± 0.010 


Extr. 


4.74 ± 0.15 


0.01 ± 0.02 
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